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Abstract 

Single transverse-spin asymmetries have been studied intensively both in experiment and the- 
ory. Theoretically, two factorization approaches have been proposed. One is by using transverse- 
momentum-dependent factorization and the asymmetry comes from the so called Sivers function. 
Another is by using collinear factorization where the nonperturbative effect is parameterized by a 
twist-3 hadronic matrix element. However, the factorized formulas for the asymmetries in the two 
approaches are derived at hadron level formally by diagram expansion, where one works with various 
parton density matrices of hadrons. If the two factorizations hold, they should also hold at parton 
level. We examine this for Drell-Yan processes by replacing hadrons with partons. By calculating 
the asymmetry, Sivers function and the twist-3 matrix element at nontrivial leading order of as , we 
find that we can reproduce the result of the transverse-momentum-dependent factorization. But we 
can only verify the result of the collinear factorization partly. Two formally derived relations between 
Sivers function and the twist-3 matrix element are also examined with negative results. 



1. Introduction 

Single transverse-spin asymmetries have been observed in many experiments[Tl[2l[3ll31[5]. To generate 
a single transverse-spin asymmetry(SSA) it requires nonzero absorptive part of scattering amplitude and 
helicity-flip interactions. This indicates that the asymmetries are T-odd effects. The study of SSA of 
hadron scattering provides a new tool to explore hadron structure and nonperturbative properties of 
QCD. In some cases like production of a heavy quark with transverse polarization, SSA can be studied 
by using perturbative QCD directly|6j. However, in the cases studied in experiment, in which an initial 
hadron is transversely polarized and is involved in the scattering, it is not possible to use perturbative 
QCD directly. To study SSA of hadronic processes two factorization approaches have been proposed. One 
is by using transverse-momentum-dependent(TMD) factorization, where one takes transverse momenta 
of partons in hadrons into account. Another is the collinear factorization. An up-to-date review about 
studies of SSA can be found in [7]. In this work we will focus on the two approaches in Drell-Yan 
processes. 

In the approach of TMD factorization, the origin of SSA arises from a correlation between the trans- 
verse spin of the initial hadron and the transverse momentum of partons in the hadron. This correlation 
is parameterized by Sivers function [H [9]. In this approach the effect of helicity-flip interactions and the 
T-odd effect are contained in the Sivers function. The helicity-flip of a initial hadron can happen because 
of orbital angular momenta of partons. This can be seen clearly in terms of light-cone wave functions [lOj. 
Hence SSA in this approach is sensitive to orbital angular momenta of partons. The T-odd effect in 
Drell-Yan processes comes from the initial state interaction, in contrast to semi-inclusive DIS where T- 
odd effects come from final state interactions. We note here that so far TMD factorization has been 
examined carefully only for physical quantities which do not contain T-odd effectsfH] [I2l [TU [HI [15] . 
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TMD parton distributions entering the factorization for these physical quantities can be defined with 
QCD operators consistently. Intensive efforts in theory has been spent to study how to consistently 
define or interpret Sivers function as a parton distribution which is gauge invariant and contains initial- 
or final state interactions [U [l6l [l8l [19]. Through these studies the role of gauge links used to define 
Sivers function becomes clear and it shows that the Sivers function in Drell-Yan processes is related to 
that in semi- inclusive DIS. The approach of TMD factorization has a simple parton- model interpretation. 
Because of this, SSA has been studied extensively in terms of Sivers functions [20l EH |22l |23l [Ml [25] . 
These functions have been also studied with models [261 1271 128[ [29j. The approach of TMD factorization 
has the limitation that it is only applicable in certain kinematic regions, e.g., in a Drell-Yan process the 
region is where the transverse momentum q± of the lepton pair is much smaller than its invariant mass 
Q. 

In the approach of collinear factorization SSA is factorized with twist-3 matrix elements |30[ I31[ 132] . 
or called ETQS matrix elements. In this approach the twist-3 matrix elements, or the corresponding 
parton distribution functions defined with twist-3 operators of QCD, contain only the effect of helicity- 
flip interactions which is taken as nonperturbative effect. The nonzero absorptive part or T-odd effect 
is not contained in the twist-3 matrix elements. It is generated by poles of parton propagators in hard 
scattering. The twist-3 matrix elements characterize the correlation between quarks and gluons inside 
the transversely polarized hadron. Therefore, in this approach SSA is sensitive to the correlation. From 
this point of view the approach of collinear factorization seems different than the approach of TMD 
factorization. However, recent progress shows that the two approaches can be unified in the kinematic 
region of q± « Q|33]. We note here that the approach of collinear factorization is applicable for the 
whole kinematical region if is enough large. The fact that the two approaches in the region of q± << Q 
indicates that there exists a relation between Sivers function and twist-3 matrix elements. Such a relation 
has been found in ^33j. There also exists another relation between Sivers function and twist-3 matrix 
elements [m [M]. Applications of the collinear factorization for SSA can be found in [35 1 [36 l [37]. 

It should be noted that in the two approaches the factorization is derived or proposed rather formally 
in the sense that one works at hadron level by using the diagram expansion. In the expansion one usually 
divides a given diagram with hadrons into various parts. Among these various parts, one consists only 
of partons. In other parts hadrons are involved. These parts represent nonperturbative effects related 
to the hadrons and they are parameterized by various parton density matrices of hadrons. It should be 
also noted that QCD factorizations, if they are proven, are general properties of QCD Green functions. 
It means that the two factorization approaches, if they hold, they should also hold by replacing hadrons 
with partons. It is the purpose of the study presented here to show how SSA in Drell-Yan processes 
can be factorized in two approaches by replacing hadrons with partons and to examine the two relations 
between Sivers function and twist-3 matrix elements. Our study is performed at leading order of as- In 
order to generate SSA in Drell-Yan processes and nonzero q±, there must be exchange of two gluons at 
the leading order. It results in that SSA at parton level is already at order of in comparison with 
the leading order of the unpolarized cross-section which is at a^. Hence it is nontrivial to show the 
factorizations at leading order. In this work we will take Drell-Yan process as an example. We replace 
the two hadrons in the initial state with a quark q and an antiquark q. In order to have helicity-flip we 
keep the quark mass as nonzero and every quantity is calculated at leading power of m. The perturbative 
coefficients in the factorization formulas do not depend on the quark mass m. It turns out that the 
proposed TMD factorization of SSA holds at the parton level, while only a part of results of the proposed 
collinear factorization for SSA can be verified with our partonic results. We also find that the two 
relations between Sivers functions and twist-3 matrix elements do not hold in general. The two relations 
need to be modified. 
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Our paper is organized as the following: In Sect. 2 we give our notations for Drell-Yan process and the 
formulas of two factorization approaches of SSA. In Sect. 3 we present our result of Sivers function with 
a parton state. SSA of Drell-Yan process with the parton state is calculated in Sect. 4. Sect. 5 contains 
the result of twist-3 matrix element with the parton state, where we show that only a part of our result 
matches the formula of the coUinear factorization of SSA. We summarize our study in Sect. 6. 



2. Notations and Factorization Formulas 

We consider the Drell-Yan process: 

hA{PA, s) + hBiPs) ^ l*{q) +X^r+£++X, (1) 

where Ha is a spin- 1/2 hadron with the spin- vector s. We will use the light-cone coordinate system, in 
which a vector a'^ is expressed as a'* = {a'^,a'' ,a±) = ((a° -|- a^)/\/2, (a° — a^)/\/2, a^, a^) and a\ = 
(a^)^ + (a^)^. We also introduce two light-cone vectors: n'' = (0,1,0,0) and = (1,0,0,0). We take a 
light-cone coordinate system in which the momenta and the spin are : 

P^B = (PXb^PabAO), s'^ = {0,0, s±). (2) 

Ha moves in the z-direction, i.e., is the large component. The spin of Hb is not observed. The 
invariant mass of the observed lepton pair is Q'^ = q^. The relevant hadronic tensor is defined as: 

^'^^ = E/ ■^e''i-^hA{PA,s±),hB{PB)\mYQmX){X\q{x)rq{x)\hB{PB),h (3) 



X 



and the differential cross-section is determined by the hadronic tensor as: 



da Anal^Q'^ 



em^q 



dQ^d^q±dq+dq- 3SQ^ 



K<l'-Q'){^-9,.)W'^^. (4) 



We are interested in the kinematical region where q\ « Q^. The hadronic tensor at leading twist 
accuracy has the structure: 



nv al3 117(1) / aa v . av U ui' Q0 \ 717(2) 

+q±a {eTql + ^Tql) q± ■ s±wP + ■■■ (5) 

with the notation: 

5f = 5'^^ - n^r - n^'l'^, = e^^'^^/^n^j. (6) 

In the above, we only give the structures symmetric in ixv. Wj:\i = 1, 2, 3) represent T-odd effect related 
to the spin. W^''^'^ are responsible for unpolarized cross-sections. W^,^^ contributes to SSA in the region 
ql « 

dQ^d^q^dq+dq- " dQ^d^q^dq+dq- = ^SQ^^^'' ' ^ ^^"^^'^^^ V + ^^^^/^ )J " 
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We will focus on w!^^ to see if it can be factorized. 



In the kinematical region >> 9^ ~ ^QCD TMD factorization can be performed for Wrp based 
on the diagram expansion. The result at tree-level can be written as a convolution with Sivers function 
and TMD parton distribution [T31 [33]: 

w!p\zi,Z2,q±) = irr d^h±d^h±^^-2^<l±izi^h±)q{z2,k2±)d^{ki± + k2± - q±)H, (8) 

where the variables zi 2 are defined as: = ziP^ and q~ = Z2P^. H is a perturbative coefficient, i.e., 
H = 1 -\- 0{as). Beyond tree- level one has to implement a soft factor representing effects of soft-gluon 
radiation. In the above q± is Sivers function. To define it with QCD operators we introduce a gauge link 
along the direction u with = {u^ ,u~ ,0, 0): 



Lu{-oo,z) 



Pexp^—igsJ dXu-G{Xu + z)^ 



t 



(9) 



The Sivers function relevant for Drell-Yan process is defined in the limit << u [9| 113^ [T8]: 

q±{x,k^)eTsL,k^u = ] f %^e-'^-'{{PA,s^mz)Ll{-oc,z)^+Lu{-^,0)m\PA,s^ 



4 J (27r)3 

-{s± ^ -s^)}, (10) 

with z^^ = {0, z~ , z±). X is defined as k^ = xP^. Beside the renormailzation scale /i, Sivers function also 
depends on the parameter: 

e = '-^(p-)\ (u) 

The limit « u~ is to be understood that we discard all contributions in Eq.(9) which are zero 
with ^ 00. The definitions of other TMD parton distributions of a unpolarized hadron can be found 
in |13[ I14j . In the TMD factorization of SSA one takes transverse momenta of incoming partons into 
account. The T-odd effect and spin-flip effect are parameterized by Sivers function. It should be noted 
that the TMD factorization can be extended to the region » q^ » Ag^^,. 

A collinear factorization can also be performed for SSA or W^^'', where the T-odd effect comes from 
poles of partons in the hard scattering and the spin-flip effect is parameterized with the twist-3 matrix 
element which is defined as|301 [31]: 

5s f ^Iy]dy2 iy2(x2-xi)P+ -iyixiP+ uu 



TFixi,X2)s'i = ^ J 



4tt 



■ {{PA,s±\^P{yln)J+G+^iy2n)^P{0)\PA,s±) - (s^ ^ -s±)} . (12) 

In the above we have suppressed gauge links between operators. These gauge links are defined with the 
vector n and make the above definition gauge invariant. One can also view the definition as given in the 
gauge n ■ G = 0. With the twist-3 matrix element SSA or takes the following factorized form|33j: 

W^^^ qfg) Hc<^Tf, (13) 

where q is the standard parton distribution. He is a coefficient function calculated perturbatively. Its 
leading order is at a^. Details about the above result can be found in [33]. It should be noted that the 
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above collinear factorization is derived for the kinematical region with » Ag^^^, and q\» Agj^^,. 
It should also be valid for the region with Ag^^, « « Q^. In this region the two factorization 
approaches apply. It has been shown both factorizations give the same results [33] in that region. Hence 
a relation between Sivers function q± and the twist-3 matrix element Tp can be found. 

As discussed in the introduction the derivation of these factorization formulas is based on the diagram 
expansion, in which one works with hadronic states by introducing various parton density matrices of 
hadrons. If these factorization formulas hold, it should also hold if we replace hadrons with partons. It 
is the task of the subsequent sections to check these factorizations and different relations between Sivers 
function and the twist-3 matrix element by replacing hadrons with partons. 
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Figure 1: Diagrams for the contributions to Sivers function. The double lines represent the gauge link. 

3. Sivers Function with a Quark-State 

We replace in the definition of Sivers function the hadron Ka with a quark q. The quark has the 
momentum p^^ = (p"*",!?", 0, 0) with = m? . The finite quark mass m will introduce effects of spin- flip. 
In order to have T-odd effect, exchanges of virtual gluons should present in the amplitude. Also, at least 
one gluon must be in the intermediate state to generate nonzero k^. With these considerations one can 
find the possible contributions to Sivers function at leading order of Ug- These contributions are given 
by diagrams in Fig.l., where the contributions are represented as the interference of amplitudes. The 
amplitudes of Fig. la to Fig.le can have absorptive parts. The interference of these absorptive parts with 
the amplitudes represented with Fig.lg and Fig.lh will give nonzero contributions to Sivers function. 

At first look there are many diagrams contributing. However, the Sivers function q±{x, k_[_) is defined 
in the limit u~ « . In this limit one can easily show that only the interference of Fig. lb and Fig.lc 
with Fig.lg are nonzero, other interferences are zero. The absorptive parts of amplitudes can be obtained 
with the standard Cutkosky cutting rule and they can be represented with cut diagrams. In Fig. 2 we 
give the diagrams for the absorptive part of Fig. lb and Fig.lc. In determining physical cuts one should 
note that the energy fiow of each particle crossing the cut should be in the same direction. In our case, 
the gauge link represents an incoming particle with an infinity-large — component of its momentum. 
Therefore, the energy flow of the gauge links in Fig. 2 are from the right side to the left side. Keeping 
this in mind we flnd only one cut diagram for Fig. lb with a physically allowed cut. For Fig.lc there are 
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Figure 2: The absorptive part of Fig.lb and Fig.lc. The broken hne is the cut. 



two possible cuts. A cut cutting a real particle has no effect. According to the Cutkosky cutting rule 
one should replace all i's of propagators and vertices with —i in the left-upper part of the cut diagrams 
in Fig. 2. In this part of the cut diagrams the joining point of the gauge link and the quark line should 
be taken as a vertex, reflecting the fact that the particle represented by the gauge link annihilates the 
particle represented by the quark line. Hence, this joining point or vertex contributes an extra minus 
sign when evaluating the absorptive part of the cut diagrams. This extra minus sign can also be verified 
by a direct calculation of the interference of the complex conjugated Fig.lb or Fig.lc with the complex 
conjugated Fig.lg. 

The interference of Fig.lb with Fig.lg can be written as: 



1 



4 rabc^ 
'tis J 2 

-i2iT5{u 



4A^c 2 



1 f d% d^q 



(2vr)4 (27r)4 (2^)<^(^p)'^(^ -P + K)^"ik± + kg±) 



-i2iT6{{p - kg - q)"^ - m?) 



{p — kgY — rn? — ie {q + kg)^ + ie q"^ + ie 

■u{p, s±)7m(7 • {P - kg) + m)'y+{-f ■ {p - kg - q) + m)-fpT''T^T''u{p, 
• [{-kg + q)Pu>' + {-2q - k^ruP + {2kg + q) ■ ugf^P] - {s^ ^ -s^) , (14) 

where kg is the momentum of the gluon in the intermediate state and q the momentum of the gluon 
emitted from the gauge link. It is straightforward to calculate the contribution under the limit n+ 0. 
We have: 



q±{x, k±)\bg 



mat 



x{l — x) 



1 



(1 — x)'^im? 



In 



(15) 



levr^^"''^ ^' k\ + {I - x)'^m? k^ \k\ + {I - x)'^m? ) ' 
The absorptive part of Fig.lc receives contributions from two cut diagrams in Fig. 2. It is easy to 
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show that the contributions from the two cut diagrams cancel each other. Therefore the only nonzero 
contribution is from Fig. lb. The final result of Sivers function is: 

, , , mai.^^n x(l — x) 1 , / (1 — x)^m^ \ 

t,) , 8^(JV,' ' + (1 - J)-m^ g [ ki + (1 - x)^m^ j ■ ^''^ 

We note that the same diagrams in Fig.l will also contribute to Sivers function in DIS, where the gauge 
link is pointing to the future. This gauge link then represents an outgoing particle with an infinity-large 
— component of its momentum, hence the energy flow along the gauge link is reversed. This will lead to 
cut diagrams other than those given in Fig. 2. It turns out that Sivers function for DIS calculated with 
Fig.l is the same as the above, except a sign difference as expected. We also point out that there can be 
a difference in calculations of Sivers function or SSA between different ways to generate absorptive parts. 
In the collinear expansion imaginary parts of amplitudes are generated by poles of parton propagators. 
In general the generated imaginary parts can not be the absorptive parts, because these poles do not 
necessarily correspond to physical cuts, although the same results may be obtained. The difference and 
SSA in DIS will be studied in a separate publication. 



4. SSA of Drell-Yan Processes 

In this section we will calculate the part of the hadronic tensor relevant to SSA. We replace Ha with 
a quark and Hb with an antiquark. We consider SSA in the process: 

q{pi, s) + q{p2) ^ l*{q) +X^r+£++X, (17) 

where the quark q is polarized with the spin vector s. By the requirement that there is a T-odd effect and 
nonzero q±, we find that at leading order of the possible contributions to SSA are given by diagrams in 
Fig. 3. These diagrams are of the hadronic tensor and the black dot indicates the insertion of the electric 
current. We denote the momentum of the gluon in the intermediate state as kg. It will be very tedious to 
obtain full results from these diagrams. However, what we need is the leading contribution in the limit 

« . It will be useful by doing the expansion in q^/Q'^ first and then to perform the loop integral. 
A convenient way for the expansion is to analysis different regions of the loop momentum. 

For doing the expansion we note that each contribution from Fig. 3 can be written in a generic form 

/ ^^(2.)S(l4)(2.)'S'(p, ■ ^ I- ■ -1 . (18) 

where k is the momentum of the virtual gluon. In each contribution there are five propagators, their 
denominators are denoted as Di{i = 1, 2, • • • , 5). The nominator represented in the above as Tr [• • •] is a 
trace of product of 7-matrices. We scale the momentum kg± = — gj_ as at order of A. and expand each 
contribution in A. It is clearly that the leading order contribution in A comes when the denominators of 
all propagators are at order of A^. The power-counting is not affected by taking a cut to cut propagators. 
The requirement that all denominators in Fig.3a to Fig.Se are at order of A^ gives the following scaling 
of loop momenta: 

k'^ = {\\ \\ A, A), ki'g = (A:+ k-,\ A), (19) 

where one of the component k'^ or k~ is at order of A^. It depends on diagrams. However, some 
contributions can not have all denominators at order of A^, e.g., the interference terms of Fig. 3a or Fig. 3c 
with Fig.3g, the interference terms of Fig.3d or Fig.3e with Fig.3f. It is easy to find a rule to determine 
which contribution is dominant. If k~ in one of Fig.Sa to Fig.Se is at order of A^, then its interference with 
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Figure 3: Diagrams for the contributions to the hadronic tensor Wj. 



(1) 



Fig.3g is not at leading order of A. If in one of Fig. 3a to Fig.Se is at order of A^, then its interference 
with Fig.Sf is not at leading order of A. We note that the leading order of the nominators starts at 
order of 0{X) or higher. With this rule we can only have the leading contributions from the following 
interference terms: Fig. 3a or Fig. 3c with Fig.3f, Fig. 3d or Fig.Se with Fig.Sg, Fig.Sb with Fig.Sf if k~ is 
at order of A^, and Fig.Sb with Fig.Sg if k'^ is at order of A^. In these contributions all denominators 
of propagators are at order of A^. Evaluating the nominator of these contributions we find that only 
the nominator of Fig.Sb interfered with Fig.Sf and that of Fig.Sa interfered with Fig.Sf are at order of 
A, other nominators are at order of A^. Therefore the leading order contributions come only from the 
interference term of Fig.Sb or Fig.Se with Fig.Sf. We note that the leading contribution comes from the 
Glauber region of the momentum of the virtual gluon. Hence its propagator can be approximated as: 



+ ie 



-kj + ie 



(20) 



From the above analysis, one can see that the diagram Fig.Sb and Fig.Se, which can generate nonzero SSA, 
are in correspondence to the diagram Fig. lb and Fig.lc, which contribute to Sivers function, respectively. 

Again, the absorptive parts of Fig.Sb and Fig.Se can be found with the Cutkosky cutting rule. They 
can be represented by the cut diagrams in Fig. 4. For Fig.Sb there is only one physically allowed cut, 
for Fig.Se there are two. Taking the interference of Fig.Sb with Fig.Sf as an example, we have the 
contribution to the relevant hadronic tensor as: 



bf 



9^ 



2m 



-i2TT6{{pi - k- kgf - m2) -i2^5{{p2 + kf 



kn 



k'^ + ie 
1 



+ ie 



Tr 



{k + kg) 
757 ■ s±{-i ■pi + m)-fah ■ (Pi 



{Pi 



h)' 



m 



+ ie 



kg) + m)Y{l ■ P2 - m)jf3{-j ■ {p2 + k) + m)-f^' 



{l-{pi-k- kg) + m)7p] • {-kg + k)Pg^f + {-2k - kg)^g^P + {2kg + k)Pg' 



} • (21) 
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Fig.3b 




Figure 4: The cut diagrams of Fig.3b and Fig. 3c for their absorptive parts. 

It should noted that one should take complex conjugation of the right part of a cut diagram in evaluating 
the absorptive part. With this in mind the black-dotted vertex in the cut diagram of Fig. 3b contributes 
an extra minus sign. Now it is straightforward to perform the expansion in A and to pick up the leading 
contribution. We have then: 



d?ki k\n 1 

(22) 



16p+P2 J (27r)2 kl + ie [k^ + kg^f + (1 - x)'- 



For the physical process, the integration over kj_ is bounded from the above because the energy-momentum 
conservation. Since we only work at the leading order of qj_, we can integrate k'^ from to oo. For 5{kg) 
we can write it with the variables x and y as: 



6{kl) = 6{2ptp^{l -x){l-y)- ql) - ^ f i ./^ . +••-, q-=yP2, q+ = xp+ . (23) 

2pjp2 (1 - x) 



where the terms represented with • • • will not contribute. Then we have: 

(1) ^ as Nl - 1 x(l - x)5{l - y) (1 - x) W 

\hf 16ql{ql + {l-xfm^) ql + {l-xfm^' ^ ' 

There are two cut diagrams for the absorptive part of Fig. 3c. At leading order of A they contributions 
cancel each other. Hence we have the total W^^^ 

2^ aiN"^-! x{l-x)6{l-y) , {l-x)'^m'^ 
W^^\x,y,ql) = m-^^ \ ^ ^2 2n 2 ^2 2 - 25 



With the assignment of power of A in Eq.(19) we find that the leading order of Wj. is of A^^. With the 
above discussion about denominators of propagators, it seems also possible that the leading contribution 
comes from the region of the loop momentum with k± « q±, i.e., the region with the assignment of 
power of A: 

k^ = {X^X\\^X'), k^ = {k+,k;,X,X). (26) 

Performing the power counting and evaluating the nominator in Eq.(18) for each diagram we find the 
leading contributions from this region to wi^^ are at order of A ^. Calculating these contributions one 

finds W^'^ = 0. This indicates that nonzero contribution from this region to starts at order of A~^. 
One can also perform similar analysis for other regions of the loop momentum. The conclusion is that 
the leading term of in q± comes from the region specified with Eq.(19). 
The tree-level TMD parton distribution is: 

q{x,k±) =6il-x)6\k^)+0{as). (27) 

Our results in Eq. (I16|25|27|l shows that the TMD factorization for Wj^^ in Eq.([8]) is verified at lead- 
ing order. Although all calculations performed here are at leading order, they are non-trivial. With 
the factorization at leading order the nonzero transverse momentum q± are generated by the nonzero 
transverse momenta of incoming partons which are q and q. Beyond the leading order, the transverse 
momentum qj_ can also be generated by soft gluon radiation, a soft factor should be included in Eq.([8]). 
This soft factor can be identified by extending our calculation to the next-to-leading order as shown in 
TMD factorization of T-even quantities [III HI ISl E] . 




(a) (6) 



Figure 5: The diagrams for the twist-3 matrix element in the n ■ G = gauge. 



5. Twist-3 Matrix Element with a Quark-State 

Now we turn to the twist-3 matrix element. With the single parton state as used for calculating Sivers 
function we can also calculate the twist-3 matrix element Tp(xi,X2) defined in Eq. (jl2p . The calculation 
can be done easily in the gauge n ■ G = 0. In this gauge the gauge links in Eq. (jl2p become a unit matrix 
in the color-space. The leading order contribution comes from diagrams given in Fig. 5. A straightforward 
calculation gives: 

f (fk± 1 

Tf{xi,X2) = 2TTasGFm{x2- xi)'^6{l - X2) J ^ _ ^ )2 ~^ ^ ^2)- (28) 
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The results are U.V. divergent. We can regularize the U.V. divergence with the dimensional regularization 
and derive the renormalized Tf(xi, X2, /x): 

2 

rp(xi,X2,/x) = ^CFm{x2 - xif5{l - X2)ln- ^ „ „ + {xi ^ X2). (29) 

2 (1 — xij^m^ 

We note that Tp(xi,X2) is zero with xi = X2 at leading order and the renormalization scale acts 
effectively cutoff of the transverse momentum. 

In the collinear factorization of SSA in Eq. ijlSh . the leading order of the perturbative function He is 
at as and He does not contain explicit /x-dependence. With the calculated Tp and leading order results 
of the standard parton distribution one can find that the collinear factorization of SSA with the twist-3 
matrix element fails to reproduce the partonic SSA in Eq. (j25p completely. In Eq. ijlSp the contributions 
can be divided into the contributions from soft-poles of propagators and hard-poles of propagators. The 
contribution of soft poles is proportional to Tp(x,x). Because Tp{x,x) with our partonic state is zero at 
leading order of Og, we can not verify the soft-pole contribution with our results. With our results one 
can actually derive a factorized formula with x < 1 and fJ- = q± as: 

Wi^\x,y,q±) = - 2(\\2 /— — g(^2)'^(l -6)7: — '^—^xTF{x,yi,qi_), (30) 

with = X / yi, S,2 = y / y2 and q{y2) as the standard antiquark distribution. There is a certain ambiguity 
to determine the hard kernel in the above. Because TF{x,x,q±) = one can replace the above -|— 
distribution with 1/(1— ^1). The above factorized contribution only corresponds to a part of the hard-pole 
contribution in the approach of the collinear factorization. The above factorization has a clear meaning 
in physics. In Fig. 3b the gluon emitted by the quark is collinear to the quark in the limit qj_ « Q^. 
This collinear system is factorized into the twist-3 matrix element Tp. Our study also indicates that the 
relation between q±{x, k±) and Tp{xi, X2) for kj_» Aqq^, derived formally in [33], is not satisfied with 
our partonic results. It will be interesting for a further study to show how one can find all contributions 
in the approach of the collinear factorization with some partonic states. 
There is another relation between q± and Tp derived formally in [19\ 



Tp{x,x)= (fk_Lklq±{x,k±). (31) 



With our result of Sivers function in Eq. ()16l) one can realize that the integral of /cj_ is U.V. divergent. 
Therefore the integral without any subtraction is not well-defined. A possible way to modify the integral 
meaningful in order to relate Sivers function with Tp is to define: 

Q±{x,b) = / d^k±k'j_q±{x,k±) exp{—ib ■ k±) 



167r 



(N^ - l)x{l - x) In^ 



(1 - x)2m262e2T 



+ 0(6^), (32) 

where we have used our result in Eq.([T6]) to produce the expression in the second line. One may expect 
that Tp{x,x) is somehow related to Sivers function in the impact space. We examine this in the below. 

We have seen at leading order Tp(x,x) = 0. This function can be nonzero at higher orders where a 
gluon can be in the intermediate state, i.e., a gluon crosses the cut in Fig. 5. These diagrams are given in 
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+ H.C 



Figure 6: Diagrams for the contributions to Tf{x,x) 



Fig. 6. For xi = X2 one can easily find that only the contribution from Fig. 6c is not zero. The calculation 
is similar as before. We find 



1 6c 



-^9si-2-x^)p^ K J WfW)"^ ^^^^" 

■■k5{xip'^ — {p^ — k't — k~^))6{k'^ — (x2 — xi)p~^)e'^ 



1 



1 



1 



1 



fc^ + ie {k + kgY + ie {p — kg) — m? — ie [p — kg — kY — m'^ + i£ 



8imk^^{-l)efk^as±p'i^ {p+ - k+ + ;-k+ ] [1 + 0{k+)] . 



k1. 



■k+ 



(33) 



Because we will encounter U.V. divergences we have used the naive 75-prescription in d = 4 — e dimen- 
sion. The term with e in the last line comes from the trace of 7-matrices. Performing integrations and 
subtracting U.V. divergence we have for Tp(xi,X2) with xi = X2 = x: 



Tf{x,x) 



maz. 



{Nl - l){l-x)x In' 



+ ln- 



(1 - xYm? ' (1 - xYm? 12 2 



{Nl-l){l-xf 



1 



+ ln' 



(34) 



47r ^ '''' \2 ' ~" (\ - xfrrfi j ' 

Prom the above results it is clear that one can in general not write down a factorized relation like 

TF(x,x,/i) = C(x,At,6) Jd^k±k'iq±{x,k_L)ex.p{-ib-k±) + 0{b) (35) 

with the coefficient C{x,fj,,b) as a perturbatively calculable coefficient. 



6. Summary 

SSA is a T-odd effect and it requires helicity-flip interactions. Two factorization approaches to study 
SSA has been suggested. One approach is by using TMD factorization which takes transverse momenta 
of partons into account. In this approach the T-odd effect and the effect of helicity-flip interactions 
are parameterized by Sivers function, which can consistently be defined with QCD operators. Another 
approach is to use standard collinear factorization. In this approach the effect of helicity-flip interactions 
is parameterized by the twist-3 matrix element, while the T-odd effect arises from hard scattering of 
partons. Sofar, all factorized formulas of SSA in the two approaches have been derived in a rather formal 
way in the following sense: One has used the diagram expansion with hadrons and has divided a given 
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diagram with hadrons into various parts. Among these various parts, one consists only of partons. In 
other parts hadrons arc involved. These parts represent nonperturbative effects related to the hadrons 
and they are parameterized by various parton density matrices of hadrons. Since a proven factorization 
is in general a property of QCD, the above two factorizations should hold by replacing hadrons with 
partons, if the two can be proven. The factorization approaches have been not examined with partonic 
states. 

The study of our work presented here is to examine the two factorizations of SSA in Drell-Yan 
processes with partonic states in first time. We replace the initial hadrons with a quark and an antiquark, 
where the quark is transversely polarized. With the quark state we can calculate Sivers function and 
twist-3 matrix clement. SSA can be calculated with the quark-antiquark state. The finite quark mass is 
introduced to flip helicities. With our partonic results we can examine the two factorization approaches 
at leading but nontrivial order of ag. It turns out that our partonic result of SSA can be factorized within 
the approach of TMD factorization, but the results of the proposed factorization formula of the collinear 
factorization can not be verified completely with our partonic results. This is our main result. Using 
our results we have also examined two formally derived relations between Sivers function and the twist-3 
matrix element. Our partonic results do not satisfy these two relations. It will be interesting to see how 
one can derive the formally derived results in the collinear factorization with some partonic results. 

In our work we have taken a quark-antiquark state to replace the initial hadrons in Drell-Yan processes. 
We can extend our study to the case of a quark-gluon state to examine the two approaches. It is worth 
to point out that corrections of the perturbative coefficient at higher orders of can be studied after the 
verification of a factorization with partonic states at leading order. For SSA the verification at leading 
order is nontrivial as shown in this work. With our progress to understand SSA in Drell-Yan processes, 
we are able to study SSA in semi-inclusive DIS and other possible processes. Such a study is currently 
under the way. 
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